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Abstract 

O . 

, A new nonperturbative approach is used to investigate the Gross-Neveu model of 

four fermion interaction in the space-time dimensions 2, 3 and 4, the number N of 
inner degrees of freedom being a fixed integer. The spontaneous symmetry breaking is 
shown to exist in D = 2, 3 and the running coupling constant is calculated. The four 
dimensional theory seems to be trivial. 

o 
o 

O 1 Introduction 

r- . 

Q\ , In spite of its great and numerous successes the perturbation theory cannot describe a wide 
class of important phenomena (like confinement or spontaneous symmetry breaking) playing 
the key role in the problem of full and comprehensive description of physical reality. So it is 
quite natural that considerable efforts are applied in order to develop the nonperturbative 
methods in the quantum field theory. So far, only a few sufficiently effective methods like 
the effective potential, the 1/N expansion, the Gauss effective potential method |l| or the 
^ \ variational perturbation theory [fj are known. 

A new nonperturbative approach has been recently proposed in work ||. The ability 
of the method was demonstrated by the example of a self interacting scalar field in various 
dimensions. In the present paper we would like to investigate the Gross-Neveu model of 
the fermion fields with an arbitrary fixed number N of inner degrees of freedom. The cases 
D = 2,3,4 were elaborated and the spontaneous symmetry breaking was found to exist in 
two and three dimensions. For the four dimensional Gross-Neveu model our consideration 
gives the arguments in favor of the triviality of this model. These results exhibit the efficiency 
of the method and they are the finite N generalization of the known results obtained in the 
framework of 1/N expansion 0-0. 

The paper is organized as follows. Section 2 consists of a short introduction to the 
method of |5| to be used throughout the paper. Sections 3, 4 and 5 are devoted to the 
Gross-Neveu model in the two, three and four dimensions, respectively. Section 6 contains 
another approach to the subject based on the bilocal source formalism. 
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2 The method 



In our approach we will use one of the most suitable tools for the nonperturbative treatment 
of a quantum field model — the Schwinger- Dyson equation. As the system of Schwinger- 
Dyson equations for the Green functions consists of the infinite number of mutually connected 
equations, one should truncate it in some way in order to find an approximate solution. It 
is obvious enough that the concrete way of truncation has a crucial significance for the 
results. For example, if we solve the Schwinger-Dyson system iteratively by expanding the 
Green functions into the series in the coupling constant, we obtain the perturbative solution. 
Being the simplest from the practical point of view the perturbative approach is the worst 
in the mathematical sense. The matter is that a small parameter (the coupling constant) 
is a multiplier at the highest derivative term of the functional differential Schwinger-Dyson 
equation for the Green functions generating functional. This means that the equation is 
a singularly perturbed one. So, the perturbative procedure is valid only for the restrictive 
class of boundary conditions and cannot, in principle catch all solutions ||. 

Taking into account the reasons above, we can conclude that a proper approximation 
scheme for the nonperturbative solution of the Schwinger-Dyson equation should obey the 
following requirements: 

• It must take into account the highest derivative term of the Schwinger-Dyson equation 
already at the leading approximation. 

• It must allow one to make the renormalization procedure. 

• It should be simple enough for the practical calculations. 

An approximation scheme proposed in 0| obeys the requirements listed above. Below we 
give a short introduction to the main ideas of the method for the reader's convenience. 
Consider the theory of a self-interacting scalar field <fi(x) with the action 

5(0) =jd X {^-{d^f - -x^y (2.i) 

The generating functional of n-point Green functions can be written as follows 

oo 

G = £ G n j n , (2.2) 

n=0 

where j(x) is a field source. The n-th derivative of G at j = is the n-point Green function 
G n . 

The Schwinger-Dyson equation for the generating functional of this model reads 

^ +9 *W) +4X w^- ,3(x)G=0 - (2 - 3) 

The central idea of the iterative scheme is to consider the last term of this equation as 
a perturbation to the leading approximation. That is we take the following "equation with 
constant coefficients" as the leading approximation: 

(" 2 + d 2 )^- + 4A^ = 0. (2.4) 
dj{x) or{x) 
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Presenting the full functional G(j) as the sum 



G(j) = EG (n) (j) (2.5) 

n 

we then write for the terms of this sum the recursive chain of equations: 

^ + dT^-- + 4A^ = ij(x)G^ . (2.6) 

The solution for (|2.4j) is sought for in the form = exp (iaj) and from ( j2.4| ) one obtains 
a characteristic equation for the function a(x). For the G {n) we put G^ 1 = P n (j)G^°\ where 



P n (j) is a polynomial in j{x) with unknown coefficients to be defined from ( |2.6| ). These 
coefficients define the Green functions of the corresponding step of the iteration scheme. It 
should be noted, that at the leading approximation one defines "the vacuum" of the model, 
at the first step the connected part of the propagator enters the game, whereas for the higher 
Green functions one can define an approximant for the disconnected part only. At the next 
steps of the scheme we calculate many-particle amplitudes and corrections to the propagator 
and so on. 

This procedure is model independent and has a regular character. The last property is 
due to the generating functional is regular at j = by definition, as its derivatives at this 
point are the Green functions of the model. Therefore the perturbation theory around the 
point j = is regular recipe for solving the Schwinger-Dyson equation. 

The renormalization procedure can be easily introduced into the scheme (see || for 
details). Shortly the renormalization can be carried out in the following few steps. 

(i) All the necessary counter terms are expanded in a sum, analogous to (|27 



[i 2 — > fj 2 + <5/i 2 0) + + . . . A — > A + 5X(p) + SX(X) + . . . etc., (2.7) 

where the subscript of a counter term stands for the approximation scheme step at which 
this counter term should be taken into account. 

(ii) Equations (|2.4|) and ( |2.6| ) are modified respectively: 



if + H) + (i + 5Z T) d2 )jj^ + 4 ( A + 6X ^jf^ = ^ 



bet 1 ) 

+ + (1 + + i(x + sx m) — 



and so on. 

(iii) The counter terms entering the scheme at the n-th step can be fixed only at the 
next step. Before the renormalization of the next step has been done, the equations of the 
n-th step are nothing but some relations among the counter terms. 

Now we are ready to go to the subject of our paper — the Gross-Neveu model of the four 
fermion interaction [Hi. 
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3 The general consideration and D=2 case 



In D. Gross and A. Neveu have investigated the model of the iV-coloured spinor fields 
with scalar-scalar four fermion interaction at large N. Most of the works devoted to the 
Gross- Neveu model also deals with the 1/N expansion ||-@. Our goal is to investigate the 
Gross-Neveu model when the number of colours iV is an arbitrary fixed integer. 
The model is defined in the Z)-dimensional Minkowski space-time by the action: 

S = J d D xU(i^-m)^+^ (V^) 2 ^) , (3.1) 

where ipk{x) is a spinor field with N isotopic degrees of freedom. The summation over the 
isotopic indices k is implicit in ( |3.1| ). The mass dimension of the coupling constant A is 
zero in the two dimensional space-time, therefore model (|3.1|) is renormalizable in the two 
dimensions. 

Let us transform our model to that of the spinor and scalar fields coupled by Yukawa 
interaction [Q: 

S cS = J d D xU{i^-m)ijj- ^0(/i 2 + <9 2 )0 + 0(#)) . (3.2) 

Model ( |3.2|) is equivalent to ( |3.1|) if one identifies (3 = /ivA, — > 0/// and takes the limit 
[L — > oo. To be more exact, we define the renormalized Gross-Neveu model as /i — > oo limit 
of the renormalized effective model ( |3.2j) , \x being the renormalized mass of (see below). 
Below we will refer to this identification as the Gross-Neveu limit. The simplest way to verify 
the mentioned equivalence is to consider the path integral representation for the generating 
functional: 

GeffO') = / V<\>Vi)Vi) exp (iS eS + i J j<fi) . (3.3) 

Now we perform the Gauss integration over the spinor field in ( |3.3|) and then apply the 
just described approximation scheme to the resulting functional G(j). However there is one 
subtle point about the case N = 1 in two dimensions. As is known, the two-dimensional 
Gross-Neveu model at N = 1 is equivalent to the Thirring model. It is known, that the 
Thirring model does not reveal the spontaneous chiral symmetry breaking, so the same 
should be true for the Gross-Neveu model at N = 1. The matter is that in the case of two 
dimensions and N — 1 (and only in this case) the character of the symmetry of the Gross- 
Neveu model changes. At the chiral limit m —>■ this model is invariant under the discrete 
transformation ip — > ^^ip. But in the exclusive case N — 1 the symmetry becomes continuous: 
ip — > e tai5 "0 . Such an effect can be easily seen by passing to the spinor components: we find 
that non-invariant terms are equal to zero as a consequence of the grassmanian nature of 
the spinor field. So we have the continuous symmetry in the two dimensions, which cannot 
be broken spontaneously due to the Mermin- Wagner- Coleman theorem. 

In order to take into account this fact explicitly one could (following, for example, ||) 
integrate over only (iV — 1) components of the spinor field in the generating functional ( |3.3| ). 
The consideration will be slightly more complicated, but no new effect will arise except 
for the multiplier (N — 1) instead of N in front of all the quantities responsible for the 
chiral symmetry breaking. We choose the more transparent way of integrating over all the 
components of the spinor field but should always keep in mind that our consideration is valid 
only for iV > 2. 
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Upon integrating over the spinor field in ( |3.3|) we find the following Schwinger- Dyson 
equation for the generating functional G(j): 



where A(x) is the free fermion propagator: 



(x,y) A (y - x) 



G - ij{x)G = , (3.4) 



(ift — m)A(x) = — 5(x) , A(x) 



dPp p 1 + m 



(2n) D m 2 -p 2 ' 

and D stands for the space-time dimension. 

Writing the generating functional as the sum G = G^ +G^ + . . ., we put for the leading 
approximation G^: 



(^ + d 2 )j— + Npjd D yTr (l - (3 A 



8j(x 



— 1 (x, y)A(y- x) 



iSj 



G^ = 0. 



(3.5) 



The solution to this equation is sought for in the form G^ = exp(ia* j), where * means the 
space-time integration and a must be constant by virtue of the Poincare invariance. There- 
fore the Green functions of the leading approximation are constants and fully disconnected. 
Nontrivial connected parts of the Green functions arise at the further steps of iterations (see 
discussion in Section 2). 
Now we have from (I3 - 



ifj?a + N/3TrR(0) = 0, 
where we adopted the convenient notation R(x): 

R(x) = J d D y(l - p*A)~\x - y) A (y) ee (A" 1 - (3a) 



The Fourier transformation of the function R reads: 

d D P „- ipx _P'+M 



R(x) 



(2tt) d M 2 - p 2 
The first step equation is of the form: 



M ee m - pa . 



(3.6) 



SG& 



Sj(x) 



Jd D yTr (l-/3A~) \x,y)A{y-x) 



i8j> 



= (3.7) 

and we substitute into this equation the following expression for the G^: 

G« = (-j *£>* j +ir * j)G (0) , 

where the function T>(x) and the constant r are to be determined from equation 
Omitting some straightforward calculations we write down the resulting equation: 
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ij(x) = (/i 2 + d 2 ) [{V * j){x) + it) - iN(3 2 J d D y{V * j)(y)Tr(R(x - y)R{y - : 
+N(3 r^-yTri?(0) - iV/3 3 J d%d D y 2 V( yi - y 2 )Tr (R(x - y x )R{y x - y 2 )R(y 2 



x 



From the above relation we immediately find the equation for r, which defines the first 
correction to a, and the equation for T>(x), which is of the main interest for us: 

(/i 2 + d 2 )V(x) - iN(3 2 J d D yV(x - y)Tr(R{-y)R{y)) = i5(x) . (3.8) 



Equations ( |3.6| ) and ( |3.8| ) are still formal due to the ultraviolet divergencies in Tri?(0) 
and Tt(R- R) and we should renormalize our model in order to give them a definite meaning. 

Let us consider now the two-dimensional case that is put D = 2. Analysis of the diver- 
gencies of theory ( |3.2j ) shows that it is sufficient to introduce only one counter term Sfi 2 ^ 
to cancel the divergencies in the equations for a and T>. Then the leading approximation 
equation (|3.6|) is modified as follows: 

Nf3 M ( A 2 \ 
a - V + H) ^ v M 2 ) ' 

where we substitute the value of the Tri?(0) calculated with a momentum cutoff A. 

To verify the possibility of the spontaneous chiral symmetry breaking we should study 
the above equation at the chiral limit m — > (M — > —(3o): 

a= ^Thk^ ]Ii ( 1+ w^)- (3 - 9) 

The obvious solution to this equation a = leads (at least up to the few first steps of 
the scheme) to the usual perturbative expansion of the generating functional with massless 
fermions. Such a solution is unsatisfactory one from the physical point of view since it 
contains the tachion states which means that the perturbation is carried over the unstable 
vacuum J3J. That is why we will concentrate on a possible non zero solution, for which we 
have: 



/i 2 + H) = ^ ln ( 1 + ^) • ( 3 - 10 ) 



Now from equation (|3.8| ) we get the following expression for the Fourier image V(p 2 ): 

V(v 2 ) = - 

{P J /i 2 + ^ 2 0) +iV/3 2 S(p 2 )-p 2 ' 

where the bare mass operator S reads: 

S(p 2 ) = -iJ TrR(p + k)R(k) . (3.11) 

Let us renormalize T>(p 2 ) in the Euclidean momentum region by the condition: 



at p 2 = -uj 2 V(-uj 2 ) = 4t • (3.12) 
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With the renormalization prescription (|3.12|) we can find the sum (/i 2 + 5ji 2 ^) and write the 



renormalized function T> and the equation for the nonzero value of a at the chiral limit as 
follows: 

iP-V) = , 2 - (, 2 + «") + ^ (/(-^) - / (£)) 0.13) 

7V/3 2 AT/3 2 / u; 2 \ 2 2 

^r + ^r / l^J =/i - w ' (3 - 14) 

where the function f{6) reads: 



f(6) = f ln{l + 6x{l-x))dx. 



o 



In the region p 2 > 0, T>(p 2 ) defines the s-channel amplitude of the two fermion scattering, 
therefore, as directly follows from ( 3.13|) , the point p 2 = 4(3 2 a 2 is a two particle threshold 



and a fermion with non zero mass m F = (3a exists in our theory. 

Now we take the Gross-Neveu limit: (3 = fj,\/X, /J, — ► oo. It can be shown that at this 
limit the constant a is proportional to the fermion condensate : 



a — ► y/X^t^ ^ (3a — ► A<0V», (3.15) 



therefore 



m F = X <ipip> . (3.16) 

On the other hand the product (ij3) 2 T>(p 2 ) in the Euclidean region p 2 = —q 2 < tends 
to the running coupling constant of the Gross-Neveu model: 



Calculating the limit explicitly we find the running coupling for the Gross-Neveu model: 



X r(h^) ~ " I ATA/„o, ^ 77, ,2/2^ • (3 ' 17) 



l + ^(f(q 2 /m 2 F )-f(u 2 /m 2 F ) 



It is obvious from ( [3.17|) that A = A r (l; ^f-). Then, denoting the massless combination 
m 2 F /uj 2 as £ we find from ( |3.14p : 



A~ x (l; = — Jl + 4£arcsinh^= , (3.18) 



7T v ^ 

which gives the final expression for the coupling constant ( |3.17| ): 



2 

xJ-^y) = , T — . (3.19) 

iV^T^arcsinh^ 

No free parameters are available except for the fermion mass mp- It is in the full agree- 
ment with the number of parameters of the initial Gross-Neveu model: one dimensionless 
parameter A is changed for one mass parameter mj? 1. 
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In the deep Euclidean region q 2 — > oo, the running coupling constant asymptotically 
vanishes 




(3.20) 



which reflects the known property of the asymptotical freedom of the Gross-Neveu model in 
the two dimensions. 



4 D=3 

The three dimensional case is interesting in the two aspects. Firstly, the four fermion in- 
teraction model (|3.1|) is not perturbatively renormalizable, since the mass dimension of the 
coupling A is negative. But the effective Yukawa theory ( |3.2| ) is still renormalizable (even 
super renormalizable) model. Therefore, our approach is a method for handling the non- 
renormalizable theory. And secondly, the dimensional regularization plays a distinguished 
role in the three dimensions, due to a typical singularity of a one-loop term of the theory 
( |3.2j ) has the form r(l — D/2) or T(2 — D/2) and, therefore, the divergencies are absent in 
this regularization if D = 3. The theory is finite and no renormalization is needed. That is 
why we make the calculations in two regularizations: the dimensional and momentum cutoff 
ones and compare results. 

The form of the main equations ( |3.6|) and fl3.8|) does not depend on the concrete value of 
the space-time dimension, so we start our consideration directly with those equations. 

Dimensional regularization 

On calculating the quantity Tri?(0) in the D = 3 — 2e dimensions we obtain from ( |3.6|) 
the equation 

which transforms at the limit m — * to the following result: 

-o. (4.D 

We have two solutions again: the trivial a = and the nontrivial a = ^Ji 3 ones. The 
nontrivial solution gives the nonzero condensate at the Gross-Neveu limit: 

Pa^\<W>=—. (4.2) 



From equation ( |3.8j ) we find 



tV- 1 (p 2 ) =fi 2 -p 2 + ( l alxJ(3 2 a 2 -p 2 x(l-x) , (4.3) 

2n Jo v 

from which in turn we get the mass of a fermion: 

47T 

m F — (3o~ = A <^>= — . (4.4) 
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In this case the running coupling in the Euclidean region p = —q < reads: 



47T \ v^pv 4m|y y 4m| + g : 
Its asymptotics in the deep Euclidean region q 2 — > oo is 

i. e., the Gross-Neveu model is asymptotically free in the three dimensions. 

The momentum cutoff regularization 

The one-fermion loop Tr(R ■ R) and the fermion tadpole Tri?(0) possess the linear diver- 
gence, when calculated at the momentum cutoff A, therefore we introduce a counter term 
5/i( ) t° cancel the divergence at the leading approximation. For the nonzero value of a we 
have the following equation 

2 , r 2 3iV/3 2 ( k , 7T 



^ + H) = -2^- [ A + • (4 - 6) 

Equation (|3.8|) gives: 



i2T V) = /i 2 + 6tf Q) -p 2 + ^ Qf 1 dxy/lPcr*-p>x(l-x) - . 

For the renormalization of the theory we take the same scheme (|3.12| ) as we did in the 
previous section, i. e., we put 



V(-u 



2\ 
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at some Euclidean momentum p 2 = —uj 2 . We note, that due to the leading approximation 
equation ( |4.6| ) the renormalized function T> does not depend on the concrete regularization 
scheme and as a consequence so does the running coupling. 

Using the renormalization prescription we can define the sum [i 2 + SfiLs and get the 
renormalized equation for the nonzero constant a from ( [4.6| ). At the Gross-Neveu limit the 
equation can be written in the form : 

sin 2 z + b(z - sin2z) = , b = —\uj, (4.7) 

on 

where we have denoted 

OJ 

z = arcsin— f z G [0, 7T/2J , 

\Ju) 2 + 4A 2 <ipip> 

and have substituted [3a — > A <ipip>. One can show that equation ( J4.7|) always has a 
nonzero solution z G [0, 7r/2]. Under the condition w < A <ij)ip> the solution is of the form 

A <ipi!>-- 



3N\ 

The analytic structure of the renormalized function T> allows us to conclude that there 
exists a fermion with the mass mp = (3a in our model and to get a running coupling constant 
in the momentum cutoff regularization: 
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3NrriF I m F ( q 2 \ . q 2 



arcsiiii 



An \ ^fq 1 V Am 2 F J y Am 2 F + q 2 



which is identical to the result of the dimensional regularization ( |4.5| ). 

The only parameter of the model is the fermion mass, which defines the strength of 
interaction at the classical limit q — > 0: 

The similar results were obtained in the framework of 1/N expansion [[/J. 



5 D=4 

The case of the four dimensional space-time is more involved due to the effective theory 
contains the divergencies in the four point Green function of the scalar field. Therefore we 
have to introduce into the action the corresponding self- interaction of the scalar field ~ 4 . 
Thus we start from the action 



S cS = J ^^(z^-m)^-^^^ . (5.1) 

Here the constants Z a represent the boson field renormalization whereas the fermion renor- 
malization multipliers are considered as being absorbed into the norm of the fermion fields 
and into the constants of the action. 

The question of equivalence of the model ( |5.1| ) (also referred to as "extended Gross-Neveu 
model") to the Gross Neveu model ( p.l|) becomes nontrivial in the case D = 4. It is evident 
that such an equivalence can take place in the nonperturbative sense only. This equivalence 
(in chiral limit) was motivated in the framework of 1 /iV-expansion (see, for example ||) 
and is based on the fact that in this limit the terms <frd 2 <j) and <p A of effective action (|5.1|) 
are irrelevant when one defines physical quantities in the critical infrared region. As will be 
seen below the same arguing can be applied in our approach too. Therefore we define the 
Gross-Neveu model in D = 4 as the Gross-Neveu limit of the renormalized model (|5.1|) 

We will use the dimensional regularization in this section. As the leading approximation 
to the Schwinger- Dyson equation of the model ( |5.1| ) we have the following: 



5j(x) 3! 6f(x) 



N{3jd\Tr (l-(3A~y\x,y)A(y 



— x 



G (0) = (5.2) 



The solution to this equation is sought for in the form = exp(zcr * j), which gives the 
connection of the leading counter terms: 



7 (o) 



(5.3) 
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The principal difference with the previous cases D = 2,3 is that we have two counter 
terms in ( p.3| ) and upon their fixing equation ( |5.3[ ) turns into an identity and does not define 
any specific value of a. We will seek for the first step approximant in the same form as in 
the previous sections and find the following result for the Fourier image T>(p 2 ): 

™-\v 2 ) = 4V + Zi 0) ^f + N^{P 2 ) - Zfv 2 , (5-4) 

where the mass operator £ is one-fermion loop ( |3.11| ) . Choose the renormalization prescrip- 
tion for T> at an Euclidean point p 2 = —u> 2 : 

at p 2 ~-uu 2 iV-\p 2 ) = /i 2 + O((p 2 + 0J 2 ) 2 ). 

This prescription fixes the value of the counter terms and the sum Z^/j 2 + zf^ 1 (3 4 a 2 /2. 
Taking into account equation ( |5.3D allows one to define all the leading approximation counter 
terms. 

Let us consider the chiral limit m — > 0. First of all it is worth mentioning that equations 
( |5.3j ) and ( |5.4| ) up to notations coincide with analogous equations of 1/iV-expansion method. 
Therefore the above mentioned arguments of || about the equivalence are valid for our 
consideration too. 

We note now that there exist two different cases which do not contradict our equations: 
the first of them is o = and the second one is o ^ 0. Let us begin with the case of the 
nonzero a and take for simplicity ui = 0. The renormalized function T> has the form 

where the function $(#) stands for the integral: 

$(0) = / (1 + x (i - x )) ln(l + 6x(l - x))dx 
Jo 

and m 2 F = j3 2 a 2 . On putting as in the above sections 

{iP?V{- q 2 ) — -iP r {q 2 ) 
we obtain from ( |5 . 5|) the expression for the Yukawa running coupling f3 r (q 2 ) 

Pr(q 2 ) = 7^ v ■ (5.6) 

Since equations ( |5.3j ) are fulfilled identically, the parameter /3 r (0) is as free as the mass nip 
is. 

Now let us investigate the behaviour of the denominator of expression ( |5.6| ) in the regions 
of small and large values of q 2 . It can be easily seen that at q 2 — > this denominator tends 
to the unity while in the deep Euclidean region q 2 — ► oo its asymptotics is — q 2 \nq 2 — > — oo. 
Due to the denominator being a continuous function there exists a finite value q 2 , at which 
the running coupling constant has a pole. This indicates to the presence of a tachion state 
in our model. The only way to save the situation is to put the free parameter j3 r (0) be equal 
to zero. 

The same is true for another phase of the model, when o = 0. The only difference is that 
we must keep an arbitrary nonzero value of the renormalization point u> in order to avoid 
the infrared singularities at the chiral limit. The running coupling reads: 
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(5.7) 



1 + 



Nu) 2 f3 r {uj 2 ) 
8tt 2 



F(q 2 /u 2 ) ' 



where F(x) = x(l — lnx) — 1. The function F(x) takes its maximum at the point x = 1: 
F mai = and is negative for all other positive x. Therefore we come to the same conclusion 
as above: 

(i) if fjr- < uj 2 /3 r (uj 2 ) , two tachion states are present and the model is contradictory. 

(ii) if < u) 2 f3 r {uj 2 ) < one tachion state is present and the model is contradictory 
again. 

(iii) if f3 T (uj 2 ) = 0, the tachion is absent but the running coupling is identical zero. 
From the analysis above it follows that the renormalized effective coupling constant of 

Yukawa interaction (pipip should be zero. This means that the renormalized model (|5.1|) 
describes a free fermion and a scalar boson with selffnteraction 4 decoupled from each other 
due to trivialization of Yukawa coupling. Further, the triviality of the theory 4 in D = 4 is 
well known and was demonstrated in various nonperturbative approaches (see for example 
and references therein). It is worth pointing out that this triviality can be shown by 
means of the same method we are using in the present work ||. In the Gross-Neveu limit 
jj, — > oo the free boson disappears from the spectrum and we come to the theory of free 
fermion. 

Note, that the triviality of four-fermion interaction in D = 4 was also found in the 
framework of 1/iV-expansion and confirmed by lattice simulations |J. Our calculation is one 
more argument in favor of such a trivialization. 

6 The bilocal source 

In this section we would like to discuss another approach to the Gross-Neveu model. The 
approach is based on the bilocal source formalism, which is more natural and informative 
than that of the preceding sections. We will directly deal with the fermion degrees of freedom 
without integrating over them. This allows one to get the full information on the fermion 
dynamics. Namely it is possible to calculate the fermion Green functions like the propagator, 
the amplitude and so on. Besides, this method is technically simpler. But as a price for 
the above advantages the use of the bilocal source encounters a problem, which makes the 
results obtained in this formalism to be not very firm from the position of the mathematical 
rigour. The problem originates from the well known fact that keeping the right Bose or 
Fermi statistical properties of the Green functions is a rather nontrivial task, when the 
bilocal source is used. We postpone the discussion of this difficulty to the end of the section. 

Now let us demonstrate how the results of the previous sections can be obtained in 
the bilocal source formalism. We consider the two dimensional iV-component Gross-Neveu 
model with the action Sgn ( ft.lp and introduce a bilocal source — the function i] a p(x,y) 
depending on two space-time points and two multi indices a and (3 including the color and 
Lorenz degrees of freedom. The n-th derivative of the generating functional 



defines the 2n-th point Green function. This is an advantage of the bilocal source, since to 
find a Green function one has to calculate half as many derivatives than when a local source 
is used. The Schwinger-Dyson equation for the functional G{rj) is of the form: 
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iX l 7 — \x 1 — \ + ( i ^x-m) ai - — - + 8 af3 5(x-y)G (6.1) 

d 2 zr] orf (xz) 



5r){3y(yz) 

where the summation over the repeated indices is assumed. The leading approximation 
obeys Q6.1J) with zero right site. Substituting = exp{Tr(?7 * A)}, we find the equation 
for the leading approximation to the fermion propagator A: 

(M - i$) A (x) = 8(x) , (6.2) 
where the mass parameter M denotes the combination: 

M = m-i\Tr A (0) . (6.3) 

Since Tr A (0) is a function of M, relation (|6.3| ) is a consistency condition, which is an analog 
of the gap equation. 

To solve the next step equation 

i\ . . + (ip-m)—z — + G {1 >=7]*— — , 6.4 
or/or/ or] Or) 

we substitute the expression 



G« = (^Tr (12) F 12 * m * m + Tr A« *t])G^ , (6.5) 



2 -(12) 

where the function F^^ {xix 2 \yiy 2 ) is connected with the leading approximant for the two- 
fermion amplitude and AW is the first correction to the propagator. The equations for F 
and A* 1 ) can be easily obtained from (|6.4j ) and (|6.5|) , so we directly write down the answer 
for the amplitude omitting the explicit form of the mentioned equations: 

F^(x lX2 \y m ) = 

dz 1 dz 2 (A(x 1 -z 1 )-A(z 1 -y 1 )) a K.(z 1 - z 2 ) (A (x 2 - z 2 ) ■ A (z 2 - y 2 )) a (6.6) 

+ ... , 

where the dots stand for the disconnected part. 
The scalar kernel K is a solution of the equation 



dy(l - iXE) (x - y)K{y) = -i\5(x) (6.7) 
and the function S(x) is a fermion loop (compare with Section 3) 

E(x) = Tr A (x) A (-x) . 

Equation (|6.7|) can be easily solved in the momentum space and the Fourier image of the 
bare kernel K turns out to be: 

K(p2) = -T^m ■ < 6 - 8 > 
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Let us consider the chiral limit m — > of our model. To renormalize the amplitude at 
this stage we should introduce only one counter term 5\(p). We require the connected part 
of the amplitude at the symmetric point in the Mandelstam variables to be: 

F conn (s = t = -co 2 ) = K{-uo 2 ) = -i\ r , (6.9) 

where s and t are the Mandelstam variables. This allows us to fix the counter term 5\(p) 
and to obtain the renormalized connected amplitude 

r '> 2 > = - itW^fw (6 ' 10) 



Now we should go back to the consistency condition (|6.3| ) which takes the following form 
at the chiral limit m — > 

M = -i(X + 5A (0) ) Tr A (0) . 

Taking into account the value of 5A(o) fixed by the renormalization scheme (|6.9|) we find 
an equation defining possible values of the fermion mass in our model. Besides the trivial 
solution M — 0, there exists the nontrivial one M = mp ^ which corresponds to the 
spontaneous symmetry breaking and gives a connection among the fermion physical mass 
rrip, the renormalized coupling A r and the subtraction point u 2 (compare with (|3.18| ): 



N\\l + —f arcsinhJ^ = - . (6.11) 



UJ 



F 



With the help of ( 6.11j) we can get the running coupling constant from ( |6.10| ) 



A r (4") = , p= , (6-12) 



1,1 N\l + ^f- arcsinh 



which is identical to ( |3.19| ). 

Turning to the discussion at the beginning of this section, where is a vulnerable place of 
the above consideration? It is contained in the first term of the Schwinger-Dyson equation 
( |6.1| ). Indeed it can be easily shown that due to the Fermi statistics of the spinor fields the 
generating functional G{rj) obeys identically the relation: 

S2 ° (6.13) 



Srjp iai (yixi)5rjp 2m (y 2 x 2 ) 5^ lCf2 (yix 2 )6r]f3 2ai (y 2 xi) 
Therefore we can write the first term of (16. II) as 



8 2 G 

and the equation thus obtained will be equivalent to (|6.1| ) from the point of view of the full 
generating functional which is a strict solution of the Schwinger-Dyson equation. But for our 
approximant such a change has a crucial consequence since it does not possess the prop- 
erty (|6.13| ). And the substitution ( |6.14j ) leads to the appearance of the term A a7 (x) A 7( g (0) 
instead of A Q/3 (a;)Tr A (0) in equation ( |6.2|) . But such an equation for A (a;) does not possess 
any solution except for the perturbative one. The possible way to resolve this contradiction 
is in the use of properly symmetrized form of the first term in equation (|6.1| ). Such a new 
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equation will be still equivalent to (6.1) for the full generating functional and will respect 
the property ( |6.13| ), when the approximate solutions is used. 

It can also be shown that the important property (|6.13| ) which is responsible, in particular, 
for the crossing symmetry is restored for the Green functions successively step by step of the 
approximation scheme. For example, all the Green functions calculated from the will 
consist of disconnected parts only and will not possess the true structure demanded by the 
Fermi statistics. But taking into account the first correction not only gives the leading 
approximation to the connected part of the amplitude but also restores the true statistical 
structure for the disconnected part of the four point Green function approximant. At the 
next steps we will successively find corrections improving the connected and disconnected 
parts of the higher Green functions. The use of the symmetrized form of the Schwinger- 
Dyson equation allows one to get the properly symmetrized Green functions directly at the 
corresponding steps of the scheme. 
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